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ABSTRACT 

A hyper-accreting stellar-mass black hole has been long speculated as the best candidate of central 
engine of gamma-ray bursts (GRBs). Recent rich observations of GRBs by space missions such as 
Swift and Fermi pose new constraints on GRB central engine models. In this paper, we study the 
baryon loading processes of a GRB jet launched from a black hole central engine. We consider a 
relativistic jet powered by w-annihilation or by the Blandford-Znajek (BZ) mechanism. We consider 
baryon loading from a neutrino-driven wind from a neutrino-cooling-dominated accretion flow. For a 
magnetically dominated BZ jet, we consider neutron-drifting from the magnetic wall surrounding the 
jet and subsequent positron capture and proton-neutron inelastic collisions. The minumim baryon 
loads in both types of jet are calculated. We find that in both cases, a more luminous jet tends to be 
more baryon poor. A neutrino-driven "fireball" is typically "dirtier" than a magnetically dominated 
jet, while a magnetically dominated jet can be much cleaner. Both models have the right scaling to 
interpret the empirical L — Li SO relation discovered recently. Since some neutrino-driven jets have too 
much baryon loading as compared with the data, we suggest that at least a good fraction of GRBs 
should have a magnetically dominated central engine. 

Subject headings: accretion, accretion disks-black hole physics-magnetic field 



1. INTRODUCTION 

Recent observations of gamma-ray bursts (GRBs) with 
the space missions such as Swift and Fermi have greatly 
enriched our knowledge of this phenomenon. These 
new data place important constraints on the GRB cen- 
tral engine models. In general, a GRB central engine 
should meet the following criteria: (1) It must be able 
to power a very powerful outflow with isotropic lumi- 
nosity at least that of the gamma-ray luminosity, i.e. 
L iso ~ (10 49 - 10 53 ) erg s" 1 (e.g. Zhang & Meszaros 
2004). (2) The jet must contain a small baryon contam- 
ination, so that the it can reach a high Lorentz factor, 
typically T > 100 (e.g. Lithwick & Sari 2001). (3) The 
central engine must be intermittent to power rapid vari- 
ability as observed in many GRBs (e.g. Fishman & Mea- 
gan 1995). (4) Since a good fraction of GRBs are followed 
by erratic X-ray flares, the GRB central engine must be 
long-lived and can power delayed activities (e.g. Bur- 
rows et al. 2005; Zhang et al. 2006). (5) In some GRBs 
(e.g. GRB 080916C), the broad-band spectra show no 
evidence of quasi-thermal emission from a fireball pho- 
tosphere (Abdo et al. 2009), suggesting that at least for 
some GRBs, the central engine has to be strongly mag- 
netized (e.g. Zhang & Pe'er 2009). 

The leading model of GRB central engine is a stellar- 
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mass black hole (hereafter BH) surrounded by a neutrino- 
cooling-dominated accretion flow (hereafter NDAF) with 
an extremely high accretion rate (e.g. 0.01 — lM /s). 
There are two main energy reservoirs to provide the jet 
power: the accretion energy in the disk that is carried 
by neutrinos and anti-neutrinos, which annihilate and 
power a bipolar outflow; and the spin energy of the black 
hole which can be tapped by a magnetic field connect- 
ing the outer world through the Blandford-Znajek (1977, 
hereafter BZ) mechanism. Both models have been exten- 
sively investigated by many authors (e.g., Popham et al. 
1999; Lee et al. 2000; Li 2000; Narayan et al. 2001; 
Di Matteo et al. 2002; Kohri & Mineshige 2002; Wang 
et al. 2002; McKinney 2005; Gu et al. 2006; Chen & 
Beloborodov 2007; Janiuk et al. 2007; Lei et al. 2009). 
Some questions remain open: For example, which mech- 
anism plays a more dominant role in jet power? Do they 
dominate in different luminosity regimes? Can one dif- 
ferentiate these different mechanisms using observational 
data? 

In view of the recent observational constraints, we plan 
to systematically investigate the GRB BH central engine 
models in detail. The results are presented in two papers. 
In this first paper, we address a fundamental problem of 
GRB central engine: the baryon loading in the jet. All 
GRB prompt emission models rely on an assumed value 
of bulk Lorentz factor T. Recent broad band observations 
have led to measurement of T for a good sample of GRBs, 
and interesting correlations between T and the isotropic 
7-ray energy and luminosity, i.e. T oc oc L" s 2 with 
ai ~ cti ~ (0.25 — 0.30), have been discovered (Liang 
et al. 2010; Lii et al. 2012; cf. Ghirlanda et al. 2012). 
It is therefore interesting to investigate whether such a 
correlation roots from the fundamental physics of GRB 
baryon loading. In the companion paper, we will investi- 
gate how the two jet mechanisms confront with the data 
of prompt GRB emission and early X-ray afterglow. 
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The paper is organized as follows. In Section 2, we 
calculate the neutrino-annihilation power and minimum 
baryon loading from a neutrino driven wind from a 
hyper- accreting BH disk. The role of magnetic field is 
ignored. By considering a range of BH spin, we sim- 
ulate 2000 GRBs for their neutrino-annihilation jet en- 
ergy (which is a proxy of the isotropic 7-ray energy) and 
the maximum Lorentz factor (defined by the minimum 
baryon loading). In Section 3, we consider a strongly 
magnetized BZ jet launched from the central BH. We 
consider baryon loading into the magnetically dominated 
jet by the pickup neutron mechanism. Free neutrons can 
penetrate magnetic field lines and drift into the jet re- 
gion. Through positron capture and proton-neutron col- 
lision avalanche, the neutrons are converted to protons 
and loaded in the jet. Our results are summarized in 
Section 4 with some discussion. 

2. BARYON LOADING IN A NEUTRINO-ANNIHILATION 
POWERED JET 

The BH central engine model with a superaccreting 
disk has been studied extensively (e.g., Popham et al. 
1999; Narayan et al. 2001; Di Matteo et al. 2002; Kohri 
& Mineshige 2002; Gu et al. 2006; Chen & Beloborodov 
2007; Janiuk et al. 2007; Liu et al. 2007; Lei et al. 2009). 
In the inner region of such hyperaccretion disk, a large 
amount of energetic neutrinos are emitted, which carry 
away the viscously dissipated energy of the accreted gas. 
If the accretion rate is large enough, cooling of the disk 
should be dominated by the neutrino emission, so that 
the disk is characterized as an NDAF. 

For an NDAF with mass accretion rate Mj gn < M < 

Mt ra ,p, advection is not important. Here Mj gn and Mt rap 
are the critical accretion rates for igniting and suppress- 
ing neutrino cooling (Chen & Beloborodov 2007). If 

M < M ign , the disc temperature is not high enough to ig- 
nite neutrino emitting reactions. If M > Mt rap , the emit- 
ted neutrinos become trapped in the disk and advected 
into the black hole. For the disk with viscosity a = 0.1, 
wc find A7i gn = 0.071MQS" 1 and M trap = 9.3M0S" 1 for 
a. = 0, and M ign = O.O21M0S -1 and M trap = l.SAfoS -1 
for a, — 0.95, where a, — J,c/ (GM 2 ) is the spin param- 
eter of a Kerr BH with mass M % and angular momentum 
J.. In this case, the total neutrino power E v can be 
approximated as 



E v = eMc z 



(1 - E ms )Mc 2 



(1) 



where e is the neutrino emission efficiency, and E ms is 
the specific energy corresponding to the inner edge radius 
r ms . The expression for E ms is (Novikov & Thorne 1973; 
Wang et al. 1998), 



EL 



4\Z-R ms - 3a. 



(2) 



where R„ 



3 /r g is the the radius of the marginally 



stable orbit in terms of r g = GM,/c 2 . We have 0.06 < 
e < 0.42 for < a, < 1. The radius R ms is expressed as 
(Bardeen et al. 1972; Page & Thorne 1974), 

i? ms - 3 + Z 2 - [(3 - Z0(3 + Zx + 2Z 2 )] 1/2 , (3) 
for < a. < 1, where Z-y = 1 + (1 - a 2 ) 1/3 [(l + a.) 1/3 + 



(1 - a.) 1 / 3 ], Z 2 = {3a 2 , + Z 2 ) 1 / 2 . We have R ms = 6.0 for 
a. = 0, and R ms — 2.3 for a, = 0.9. 

The neutrino annihilation (vv — > e + e~) process can 
launch a relativistic jet reaching the GRB luminosity. 
For a system with black hole mass M, and spin a,, 
the neutrino annihilation power E„ p from the NDAF de- 
pends on the accretion rate M. For M- lgn < M < Mt rap , 
the neutrino annihilation power can be approximated as 
(Zalamea & Beloborodov 2011), 



E vj> ~ 1.1 x 10 52 



-3/2 



m 9 / 4 erg s 1 , 



(4) 

where m = M,/M Q , and m = M/MqS -1 . 

Neutrino heating in the atmosphere just above the disk 
surface results in mass-loss from the hyperaccreting disk. 
The dominant heating processes are electron neutrino 
absorption on baryons (p + v e —¥ n + e + and n + v e — > 
p + e~). For an unmagnetized neutrino-driven wind, the 

mass-loss rate M„ can be estimated as (see also Metzger 
et al. 2008; Qian & Woosley 1996) 

M v ~ lO- 6 ^/ 3 2 <e 2 o ) 5 /3 r .5/3( ro / 3) -2(/ l / r .)-iM 0S - 1 . 

where h is the half-thickness of disk, r is the disk radius, 
e u = e io x 10MeV is the energy of neutrinos, (e 2 ,) — 
13.8(kT) 2 , and T is the disk temperature (Di Matteo et 
al. 2002). Hereafter, the convention Q n = Q/IO™ in cgs 
unit is adopted if not otherwise defined. 

The disk temperature T and height h can be obtained 
by solving the set of equations describing NDAF. Since 
we are interested in the inner disk with moderate high 
accretion rate Mj gn < M < A/trap, the disk is quite dense 
and hot. As a result, cooling by pair capture on nucleons 
should dominate over v — D annihilation. In the inner 
disk, the gas pressure P gas dominates (Di Matteo et al. 
2002). The energy and angular momentum equations as 
well as the equation of state can be simplified as (e.g. 
Reynoso, Romero & Sampayo 2006; Lei et al. 2009; Liu 
et al. 2010): 
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(8) 



where p is the disk density, P is the total disk pressure, 
A nuc is the mass fraction of free nucleons. In the in- 
ner disk, it is found X nuc ~ 1 (Popham et al. 1999). 
Hydrostatic equilibrium in the vertical direction leads to 
a corrected expression for the half thickness of the disk 
(Riffert & Herold 1995; Reynoso, Romero & Sampayo 
2006). 



/ Pr 3 B 
pGMC' 



(9) 



In Equations ©-(HJ), the relativistic correction factors 
for a thin accretion disk around a Kerr BH are given by 
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Riffert & Herold (1995), i.e 

A = 1 - 2R- 1 
(7=1- 4a.iT 3/2 



a 2 i? -2 , B 



i - 3ir 



3a 2 iT 2 , 



D 



6x + 80.x 1 / 2 



3x + 2a.x 1 /2) 



2a.iT 3/2 , 
(10) 



where i? = r/r g is the disk radius in terms of r s . 
Combining Equations (JB]) - (fT0|) . we get 



T ~ 1.2 x 10 11 A°- 3 J B-°- 3 C7- 01 a - 2 TO- - 2 i?- - 3 K, (11) 



ft ~ 1.6 x lO 4 ^ - 15 ^ - 3 ^- - 5 ^ -^ - 9 ^ 1 - 3 ^ 



(12) 



A neutrino-annihilation-powered jet has an opening 
angle of vO ~ 0.1 (Aloy, Janka & Muller 2005; Harikae 
et al. 2010). Considering only the neutrino wind that 
enters this funnel, the baryon loading rate of the jet can 
be estimated as 



= 7.0 x lO- 7 ^ - 85 ^ 1 - 35 ^ - 22 ^ 2 - 



R„ 



0.32 



'"(f) 



-0.9 



0.57 1.7 
-l a -l e -l 

0.32 



M s"(l3) 



where £ = r/r ms is the disk radius in terms of r ms . 
Let us define a dimensionless "entropy" parameter 



n 



(14) 



2 : r 



where £J m = E u p + M^ v9 cr is the total matter energy 
outflow luminosity. If most neutrino annihilation energy 
is converted into kinetic energy of baryons after accelera- 
tion, and the jet would reach a Lorentz factor r max ~ 77. 

The value of the parameter 77 likely changes during a 
GRB, since the BH has spin evolution during the hy- 
peraccretion process. Without magnetic fields, the dom- 
inant mechanism is spin-up due to accretion. The process 
can be delineated by 



dM.c 2 
dt 
dJ. 

~dT 



Mc 2 E„ 



= ML, 



(15) 



(16) 



where L ms are the specific angular momentum corre- 
sponding to the inner most radius r ms of the disk, which 
is defined as (Novikov & Thorne 1973) 



Lr 



GM. 2{3VR~-2a.) 

ms 



(17) 



Since a. = J.c/ (GM 2 ), by incorporating the above two 
Equations (p~5|) and (|16)) . the evolution of the BH spin 
can be expressed by 



da, 

~~dt 



ML ms c/(GM 2 ) - 2a.Mc 2 E ms /(M m c 2 ). (18) 



Considering spin evolution, one can define an average 77 
during the evolution of a GRB, i.e. 



J E m dt 
J M- hvP c 2 dt' 



(19) 



For a hot fireball, the 77 parameter is related to the 
terminating Lorentz factor if 7/ is not too high (e.g. 
Meszaros & Rees 2000). Recent observations of GRBs 
have led to constraints of GRB Lorentz factor for a sam- 
ple of GRBs (e.g. Liang et al. 2010; Lii et al. 2012 
and references therein). By constraining F of about 20 
GRBs through modeling the deceleration bump feature 
in the early afterglow lightcurves, Liang et al. (2010) 
discovered a tight correlation between F and i? 7i i SO , i.e. 
T ~ 182(£ 7:iso /10 52 erg) - 25 Q. Lii et al. (2012) con- 
firmed the F — E 1} [ so correlation (Liang et al. 2010) with 
an extended sample (about 50 GRBs) by applying more 
methods to constrain T. They also discovered an even 
tighter correlation T ~ 249i°|° o 52 , where L lt i so is the 
mean luminosity of the burst. In Lii et al. (2012), we 
have proposed that T oc £ 7 ' 3 ,? can be explained within 
the BH-NDAF GRB central engine model using a simpli- 
fied model. Here we give much more detailed modeling 
by including the effect of BH spin. 

In FigJU we show the simulated 2000 GRBs with ran- 
dom values of BH spin a., BH mass m, accretion rate 
771 and disk mass m^. Other parameters take the typ- 
ical values. We allow a,, m, m and ma to randomly 
vary in the range of (0.1, 0.998), (3, 10), (0.01, 3) and 
(0.1, 30) respectively. We adopt a logarithmic distribu- 
tion for the accretion rate, while linear distributions for 
other parameters. 



500 




Fig. 1. — rjvs. the central engine output power E vv . We simulate 
2000 GRBs with random values of BH spin, BH mass, accretion 
rate and disk mass. The best fit gives fj oc E®-^ 7 . 

By fitting FigfTJ we find a best-fit correlation fj oc 
E®- 27 . Since 77 values are typically lower than several 
100s, it is reasonable to believe that 77 is essentially the 
bulk Lorentz factor V, so that V ~ fj oc E®- 21 . This result 
is nearly the same as what we got in Lii et al. (2012), in 
which the index is 7/27. 

7 Ghirlanda et al. (2012) confirmed the positive correlation be- 
tween r and -B 7i i so . However, since they used an abnormal method 
to estimate T, the power index of their correlation is somewhat 
larger. 
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In order to compare with the observations, we need to 
consider two effects. One is the 7-ray radiation efficiency 
77 7 . The other is the beaming effect. For simplicity, we 
assume a relatively constant rj 7 , so that L 7 = r/jE^p oc 



horizon. It can be rewritten as 



al. 2004), one may obtain a relation between and 



E ltiBO , i.e. f b oc E Jtiso 

-0.145 



0.145 



Since L 7i i SO oc E 7 ^ so , one 



can get fb oc L 7 iso , which is very insensitive to £ 7 .iso 
and E^iso- We then obtain the relation between Lorentz 
factor T and the isotropic luminosity £ 7 .iso based on the 
above scalings, i.e. 



T oc E - 27 oc (/ b L 7 , iso ) ' 27 oc L 7 ; 



0.23 
7, iso 1 



(20) 



which agrees well with the statistical correlations ob- 
tained by Liang et al. (2010) (r oc (£ 7 iso ) - 25 ) and Lii 
et al. (2012) (r oc L°f ). 

We'd like to caution that despite of the correct power 
law index of the correlation, the absolute values of rj (and 
hence, T) is typically a little bit too low. Inspecting FigQ] 
one can see that for E ut , = 10 51 erg s _1 , the Lorentz fac- 
tor T is < 40. In order to reach large F values, one needs 
to appeal to model parameters that invoke large spin and 
high accretion rate. Another possibility to reconcile with 
the data is to assume a relatively small r] 7 , so that for an 
observed L 7i i SO , the corresponding E v „ is much larger, 
which corresponds to a higher V to be consistent with 
the data. 

3. BARYON LOADING IN A BLANDFORD-ZNAJEK JET 

With magnetic field lines threading the horizon of a 
Kerr black hole, the rotational energy of the black hole 
can be extracted by the BZ mechanism (Blandford & 
Znajek 1977). The BZ jet power from a BH with mass 
M. and angular momentum J. is (Lee et al. 2000; Li 
2000; Wang et al. 2002; McKinney 2005; Lei & Zhang 
2011) 



E B = 1.7 x l0 50 alm 2 B 2 15 F(a.) erg s~\ 



(21) 



where B. A5 = B./10 15 G and 

F(a.) = [(1 + q 2 )/q 2 ][(q + 1/q) arctang - 1] (22) 



here q = a./(l + y/1 - al), and 2/3 < F(a.) < n - 2 for 
< a. < 1 . It apparently depends on M, , B t , and a. . A 
strong magnetic field of ~ 10 15 G is required to produce 
the high luminosity of a GRB. 

As the magnetic field on the BH is supported by the 
surrounding disk, there are some relations between B, 
and M . Considering the balance between the magnetic 
pressure on the horizon and the ram pressure of the in- 
nermost part of the accretion flow (e.g. Moderski et al. 
1997), we have 



n ■* ram f J ^ * 9 

<57T Airrz 



(23) 



E uP , and T oc L . The isotropic luminosity L 7 ,j SO and 
L-y are connected through the beaming factor /& <C 1, 
i.e. L-y = /tL 7 ,iso- By combining the Amati relation 
E' p oc E°-fJ Q (Amati et al. 2002, 2008; Amati 2006) and and 
the the Girlanda relation E 1 oc (E' p ) 3 / 2 (Ghirlanda et 



B. ~ 7.4 x I0 w mm- 2 {1 + y/l - aiy 2 G. (24) 

Inserting it to Equation (j2"Tj) , we obtain the magnetic 
power as a function of mass accretion rate and BH spin, 
i.e. 

E B = 9.3 x l0 53 a 2 m mX(a.) erg s" 1 , (25) 



X(a.) = F(a.)/(1 + v/l - ^) 2 - 



(26) 



It is found that X(0) = 1/6, and X(l) = vr - 2. In 
general, a faster BH spin is more favorable for GRB pro- 
duction, as revealed also by recent GRMHD numerical 
simulations (Nagataki 2009, 2011). 

We next consider baryon loading in a BZ jet. The BH 
magnetosphere has a "floor" charge density defined by 
the force- free condition (Goldrich & Julian 1969). We 
use it to define a minimum baryon loading rate of a BZ 
jet. It is given by 

M G j = 7 x 10~ 2 (m p + <;m e )4nrlB m n F c/27r 

= 2.9 x 10~ 16 (1 + 0.5<; 3 )a.m 1/2 M Q s~\ (27) 

where c/3 = c;/10 3 , and q is multiplicity of electron- 
positron pairs, which is rather uncertain. Here Op = 
0.5f2, is usually taken to maximize the BZ power, and 



n. 



a. 



gm. 2(1 + yr 



(28) 



where r. = (1 + — a%)r s is the radius of the BH 



is the angular velocity of the BH horizon. 

Since the magnetic field is supported by the accretion 
disk, in the BZ model, the BH must be surrounded by a 
hyperaccreting NDAF. The physical processes discussed 
in the previous section must still happen, which tend 
to load baryons into the jet. The main difference is that 
the magnetic field threading the BH makes a strong mag- 
netic barrier that prevents charged baryons (protons) to 
enter the jet, making a baryon poor jet (e.g. Li 2000). 
This baryon-poor jet is surrounded and collimated by an 
optically thick baryonic outflow from the hyperaccret- 
ing disk (Eichler & Levinson 1999; Levinson & Eichler 
2003). Hereafter we assume that the dominant source 
for the baryons is the neutrino-driven wind from the hy- 
peraccreting disk. Other baryon loading processes may 
happen, e.g. baryon contamination from the sideways by 
instabilities during the propagation of the jet, or baryons 
entrained from the magnetic loopes erupted from the disk 
that may enter the jet region (e.g. Yuan & Zhang 2012). 
So our model gives the minimum baryon loading in a 
magnetized BZ jet. 

A strong magnetic field may change M„ by altering the 
neutrino heating and cooling rates in the hyperaccreting 
disk (Zhang & Dai 2010). The most important effect is 
that electrons and positrons participating in the charged- 
particle reactions are restricted into discrete Landau lev- 
els (Duan & Qian 2004). For first order estimation, in 
this paper we neglect these effects. 

Since the magnetic field only affects the charged parti- 
cles, it is important to first study the composition of the 
wind. According to Pruet, Woosley & Hoffman (2003) 
and Chen & Beloborodov (2007), a hyperaccretion flow 
is neutron-rich in its inner region. The fraction of pro- 
tons f p is only around 0.1, so that the majority of baryons 
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are neutrons (Chen & Beloborodov 2007). The neutrino- 
driven wind would take the similar mass composition. 
The number density of neutrons in the wind can be there- 
fore expressed as 

Mufn 



ATrr^v w m p 



3.5 x l^A UM B- l - M C u - JJ f n a^i'€t:iR 

-0.9 32 

2 



• 1.7 ( m \ 



^ii^cm- 3 , (29) 



where f n = 1 — / p is the fraction of neutrons, e_i = e/0.1, 
/3 w ,-i = /3 W /0.1 and w w = /3 w c is the wind speed. These 
neutrons can penetrate magnetic field lines and freely fill 
any location above the disk. 

Protons are different. Because of the existence of mag- 
netic fields, only the neutrino-driven outflow in preferred 
directions, i.e. almost align with the magnetic field lines 
from the field line foot on the disk, can be ejected into 
the atmosphere. Those protons with an ejected direc- 
tion larger than an angle 9b with respect to the field 
lines would be blocked. For a rough estimate, the pro- 
ton density in the region where a local field line connects 
with the disk (which is relevant for regions outside the 
BZ jet) can be estimated as 

M u f p el 



The neutrons that enter the jet are not accelerated 
magnetically. In order to be loaded in the jet, neutrons 
should be converted to protons. The first mechanism 
would be free neutron decay. The rest-frame decay time 
scale is idccay ~ 900s, which corresponds to a typical de- 
cay radius of ~ I0 15 cm. This is not an effective mech- 
anism to load baryons near the central engine. Below 
we consider following two mechanisms that can quickly 
convert a significant fraction of neutrons to protons. 

The first mechanism is positron capture. As 
shown in Section 2, a hyperaccreting disk produces 
a strong neutrino/ anti- neutrino wind that deposit 
electron-positron pairs in the magnetized jet via neu- 
trino annihilation (vv — > e + e~). This leads to proton 
production via 

e + + n — > p + D e . (33) 

Dividing the neutrino annihilation power E vi > (Equation 
0) by the average neutrino energy (e„) ~ 10 MeV, we 
obtain a rough estimate of the number rate of e + e~ pairs: 



N, 



7 x 10 56 s -1 . 



The number density of pairs is 
around n e + e - ~ N e+e - / (irr 2 c) — 7.3 x 10 31 cm~ 3 . 

The rate of e + capture can be derived from the stan- 
dard electroweak theory (e.g., Shapiro & Teukolsky 1983; 
Bruenn 1985). In the jet region, the number density of 
nucleons are low enough to satisfy the non-degenerate 
condition. One then obtains the positron capture rate 



p 47rr 2 w w m p 
-3.5 x lO 15 ^ - 85 ^ 1 - 3 ^ - 22 /?,-^ 



0.57 1.7 
,-2«-l e -l 



n e + n = Kn n I f + (uj-Q)(w-Q) 2 , 
iQ+i 



'1- 



1 



(w - Qf 



j 2 doj, 



R 



0.32-1.7 ( m \ 



n.9 /£ 
2 



0.32 



(34) 



where u) is neutrino energy in units of m e c , Q — (m r 



r z 2 uPw!-i cm 3 ( 3 °) m p )/m e = 2.531, and K ~ 6.5 x lO^s -1 . The function 



Free protons and neutrons in the wind are coupled by 
nuclear elastic scattering. At the temperatures of in- 
terest, the corresponding rate is (er e /u) ~ 10~ 15 cm 3 s _1 , 
independent of the center-of-mass energy. For a neutron, 
the optical depth for elastic scattering with protons is 



— Q) is the Fcrmi-Dirac distribution 

1 



f+(u-Q) 



(35) 
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7"np — Tlp&elVz 

B — 2" — 1 
-0.8 °- 47 
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1.7 n0.47 
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• 1.7 ( m \ 



exp[((w -Q)- M+)/^e] + f : 

where 9 e = kT e /m e c 2 and /i+ is the positron chemical 
potential in units of m e c 2 . 
At fi + < 6 e and 6 e > Q + 1, equation p4]) is simplified 

as 

7rr 4 {2(1 -QY 
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(36) 



During the proparation of the jet, neutrons drift from 
sideways into the jet. The flux of neutrons diffusing 
into the magnetized jet is Jv{r) = \ np v n dn n /dx = 
^npV n (n n /l), where A np = l/(n p cr e i) is the mean free 
path of n — p collisions, x denotes the cylindrical ra- 
dius, I ~ (A np u„i exp ) 1//2 denotes the gradient length 
scale, t cxp — r/v n is the wind expansion time, and 
Vn = PnC = (fcT/rap) 1 / 2 is the neutron thermal speed. 
For T = 10 11 K, we have v n ~ 0.1c. For a typical jet 
opening angle 6*bz ~ 0.1, the neutron drift rate into the 
jet is 



where £(5) = 1.037 is the Riemann £- function. Here, 
we neglected the next-order terms 0(Q 2 /9 2 ), 0([i 2 /9 2 ), 
and 0[(Q + l) 5 /^] an d used the formula f£°(exp(x) + 
l)- 1 x n dx = (1 - 2") x r(n + l)C(n + 1) with T(n + 
1) = n\ for integer n. The above equation can be further 
simplified (in zero order in n/9 e ) as 

h e+n ~ 1.5 x lO^V^cirrV 1 . (37) 

The temperature of electrons produced by neutrino- 
antincutrino annihilation is given by 



1 



-aT^irr 2 /3jC 



E vv , 



(38) 



M n = 2n9Bzr 2 JD(r) 

-3.5 10 7 A°- 58 B : 



0.55 

r^Mos" 1 . (32) 



where /3j is the jet velocity. 
The timescale of neutron capture can be estimated by 

6 



Jl e + r 
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The capture time is sensitive to the distance r z from the 
BH. Close to the BH, a large fraction of neutrons that 
drifted into the BZ jet would be captured. This fractions 
drops quickly with increasing r z , since the density and 
temperature of electrons are lower at larger distances. 

The second, maybe more efficient mechanism to con- 
vert neutrons to protons is proton-neutron inelastic col- 
lision avalanche (Levinson & Eichler 2003). Protons 
entrained in the jet would be accelerated magnetically 
and soon reach an energy large enough so that inelas- 
tic collisions with neutrons would happen. These col- 
lisions (pn — > ppn + ..., pn — > ppn~) efficiently convert 
neutrons to protons (and other way round) so that the 
proton and neutron fractions become comparable. A pro- 
ton produced through positron capture would generate 
more protons via inelastic collisions with the neutrons. 
The proton fraction thus grows exponentially in what 
we term as a collision avalanche, until a proton-neutron 
equilibrium is reached. The optical depth for an inelastic 
collision of a picked-up proton with the target neutrons 
is 

t„_ p ~ a n - p r z M n / {TT0^ z r 2 z Vjm p ) 

-1.8 x 10 6 ^°- 58 B a83 /p; a i^BZ,-2 

^b, 1 ^" :? 8 ^-? 3 ^!^!)- - 55 ^?^- 1 . (40) 

where a n - p = 40 mbarn for inelastic collision (Hagiwara 
et al. 2002). 

The remaining free neutrons would decouple from pro- 
tons at a larger radius (Derishev et al. 1999; Meszaros 
& Rees 2000). The free neutrons would decay at a larger 
radius (e.g. R ~ 10 15 cm) and eventually picked up 
by the jet, and leaves some observational signatures in 
the early afterglow phase (Beloborodov 2003; Fan et al. 
2005). Eventually, all the neutrons drifted into the jet 
are loaded in the jet. In the following discussion, we es- 
timate baryon loading rate in a BZ jet as the neutron 
drifting rate into the jet, i.e. A/^bz — M n . 

For a magnetized central engine, one can define a pa- 
rameter 

_ E E m + Eb , , , > 

^ = Ta ? = - = ? Ki + ^o), (41) 

where E m = E vD + M^bzc 2 , and <r = E B /E m . 

The detailed acceleration process of a jet with both 
thermal power and magnetic power has not been stud- 
ied in detail. In general, thermal acceleration proceeds 
much faster than magnetic acceleration, so that the ini- 
tial matter power E m would be quickly converted to a 
matter flux, so that a® carries the usual definition of the 
ratio between a Poynting flux and a matter flux. If no 
magnetic dissipation occurs, the parameter 

M = ^o = ??(l + <To) = r(l + (7) (42) 

remains a constant as magnetic acceleration proceeds, 
r continues to increase while a drops (e.g. Komis- 
sarov et al. 2009; Tchekhovskoy et al. 2010). An ef- 
ficient magnetic acceleration terminates when a causal 
contact condition is broken. This occurs as the bulk 
flow Lorentz factor reaches the Alfven Lorentz factor 
T A = (1 + cr) 1/2 ~ cr 1/2 . According to Equation (|42]), 



2/3 1/3 

this corresponds to a ~ /x and T ~ /i . Combining 
thermal acceleration and efficient magnetic acceleration, 
the outflow would quickly reach a Lorentz factor 

r = max(/^ /3 ,7?) (43) 

at a radius ro ~ 2 x 10 n cm, which is not far beyond the 
radius of progenitor's envelope. 

The acceleration behavior of the jet beyond ro is sub- 
ject to uncertainties. For a non-dissipative, steady mag- 
netized outflow, the acceleration is rather inefficient (e.g. 
Lyubarsky 2010). The most efficient acceleration may be 
proceeded as 

r~r (r/r ) 1/3 , a ~ a (r/r o r 1/3 (44) 

if there is a continous magnetic dissipation due to cur- 
rent instability (Giannios & Spruit 2006) or the pulse of 
the shell is short enough to undergo an "impulsive" ac- 
celeration (Granot et al. 2011). Depends on the initial 
Ho, the jet may or may not reach the full Lorentz factor 

r max = Mo- (45) 

Rather, more likely the unsteady jet (as manifested by 
the erratic lightcurve behavior) would undergo internal 
collisions, which lead to distortion of magnetic configura- 
tions and trigger an Internal-Collision induced MAgnetic 
Reconnection and Turbulence (ICMART) avalanche to 
discharge the magnetic energy (Zhang & Yan 2011). 
Such a process would make the outflow to reach a termi- 
nating Lorentz factor T that satisfies 

r < r < r max • 

(46) 

with the explicit value depending on the detailed dissi- 
pation process. 

The central engine parameters evolve with time during 
a GRB, since the BH would be spun-up by accretion 
while spun-down by the BZ mechanism. The evolution 
equations of a Kerr BH in the BZ model can be written 
as 

-^ = Mc 2 E ms -E B , (47) 

^.=ML ms -T B . (48) 
The evolution equation of the BH spin is then 

^ = (ML m - T B )c/(GM 2 ) - 
at 

2a.(Mc 2 E ms - Eb)/{M,c 2 ), (49) 

where Tb is the total magnetic torque applied on the BH, 
i.e. 

T B = -2- = 1.7 x 10 i5 alm 3 B 2 15 F(a.) g cm 2 s~ 2 . 

SZf 

(50) 

To delineate a GRB, we average the parameters over 
time. One may define JIq = J Edtj J M h Bzc 2 dt. We can 
calculate Tq and r max as discussed above. 

In FiglH we show the simulated 2000 GRBs in the 
same way as we did in Fig. [TJ The same distributions 
for other parameters (BH spin, BH mass, accretion rate, 
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disk mass, etc) have been adopted. Both r max and Tq 
have been plotted. 




lV 10 50 10 51 10 52 10 53 10 54 
£(erg s" 1 ) 

Fig. 2. — r max = Ao (black) and Fq (gray) vs. the central engine 
output power E within the framework of the BZ jet scenario. Two 
thousand GRBs with random BH mass, BH spin, accretion rate, 
and disk mass have been simulated. Following scaling correlations 
are found: T max 0^ E 32 and To oc E 0M . The Lorentz factor of 
the GRB during the prompt emission phase is expected to satisfy 

T < r < 

1 max ■ 

By fitting Fig|2j we find the correlations r max oc £lq oc 
E 32 and T oc E 0M . In the BZ model, the relativistic 
jet dissipates its magnetic energy via ICMART with an 
efficiency ?7icmart to produce gamma-ray emission, i.e., 
L 1 ~ '7icmart-E'- The efficiency can be approximated 
as (Zhang & Yan 2011) ?7icmart ^ 1/(1 + ew), where 
(Tend is local magnetization parameter after the ICMART 
event. Typically one has ?7icmart > 50%. So the ob- 
served 7-ray luminosity is a good proxy of the jet power 
E. Again considering the beaming factor correction, one 
can derive 

r max oc £ 30 oc (/ 6 i 7 , iso ) - 30 oc L° f Q and T oc L°f Q . 
In reality, the real T should satisfy eq. (|46)) . Since both 
limiting Lorentz factors show a correlation similar to the 
observations (Liang et al. 2010; Lii et al. 2012), we sug- 
gest that the BZ model would also give a L— L correlation 
that is generally consistent with the observations. Also 
the BZ jets are much cleaner than the neutrino-driven 
ones, which overcome the difficulty of the neutrino-driven 
jets that have too much baryon contamination in the jet. 
The relatively large values of r max are not a big concern, 
since there could be other baryon loading processes be- 
sides the one considered here that would contaminate the 
jet even more. 

4. CONCLUSIONS AND DISCUSSION 

In this paper, we studied the baryon loading problem 
of a GRB jet launched by a hyper-accreting BH central 
engine. We considered two types of jet launching mech- 
anisms: the non-magnetized i/P-annihilation mechanism 
and the strongly magnetized Blandford-Znajek mecha- 
nism. For both models, we considered baryons in a 
neutrino-driven wind from a hyperaccreting disk. For 



the ^-annihilation model, the baryons launched in the 
neutrino-driven wind are mixed with the photons and 
electron-positron pairs produced by ^-annihilation, and 
thermally accelerated to reach the termination Lorentz 
factor r max ~ r\ (see definition in Equation |14|). For 
a BZ jet, on the other hand, protons are blocked by 
the strong magnetic fields at the jet boundary. Only 
a fraction of neutrons can drift into the jet. We consider 
positron capture and proton-neutron inelastic collision 
processes and argue that abour half of neutrons drifting 
into the jet can be converted to protons and be picked up 
by the jet. The other half of free neutrons would decay in 
the jet at larger radii, so that eventually all the neutrons 
can be added to the baryon loads of the jet. We calcu- 
lated the minimum baryon loading of these magnetically 
dominated BZ jets in terms of the parameter /j,q (see def- 
inition in Equation |41j). Since magnetic acceleration is 
inefficient, the final Lorentz factor of the GRB can be 
between Tq (see definition in Equation [55]) defined by 
initial efficient acceleration and T max = /io. 

A phenomenological correlation between GRB Lorentz 
factor measured in the deceleration phase and GRB 
isotropic gamma-ray luminosity has been discovered re- 
cently (Liang et al. 2010; Lii et al. 2012), i.e. V ~ L?- 30 . 
If GRB radiative efficiency does not sensitively depend on 
jet luminosity, this correlation would become a require- 
ment for any GRB central engine model. With Monte 
Carlo simulations, we have shown that both baryon- 
loading models can give cleaner jets at high luminosities. 
The slope of dependence is consistent with the observa- 
tions (see also Lii et al. 2012 for a simpler ^-annihilation 
model) . The normalizations of the correlations are quite 
different for the two models. The ^-annihilation jets are 
much dirtier than the BZ jets. For a typical jet open- 
ing angle 9j = 6 ui > ~ 0.1, we found that the resulting r\ 
in the neutrino-driven jet is typically below a few hun- 
dreds, and is only a few 10s for typical GRB luminosi- 
ties. These values are too small to be consistent with the 
GRB data. The BZ jets, on the other hand, are much 
cleaner. The r max values are typically in the 10 3 — 10 4 
range. Since the ICMART mechanism (Zhang & Yan 
2011) can efficiently dissipate magnetic energy and pre- 
vent accelerating the jets to r max , and since there could 
be additional mechanisms to load more baryons in the 
magnetically dominated jets, we argue that the magnet- 
ically dominated central engine model is more appealing 
to interpret the GRB phenomenology. In view of the low 
normalization of the T — L relation in the ^-annihilation 
model, we suggest that at least a good fraction of GRBs 
should have a magnetically dominated central engine (see 
also Meszaros & Rees 1997; Wang et al. 2002; Lei et al. 
2009; Yuan & Zhang 2012). 

In this paper, we did not consider the dependence 
of M v on the strength of magnetic fields. In strong 
magnetic fields, electrons and positrons participating in 
the charged-particle reactions are restricted into discrete 
Landau levels, which may alter the neutrino heating and 
cooling rates, and therefore change M v . We will study 
these effects in future work. 

Besides the discussed two mechanisms to launch the 
jet, it is possible that an intrinsically episodic jet is 
launched from the disk through a magnetic process (Yuan 
& Zhang 2012). The baryon loading process of this mech- 
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anism is more difficult to calculate, since baryons from 
the disk can be directly entrained in the magnetic bub- 
ble and escape. We do not discuss this mechanism in this 
paper. 

Overall, we restrict ourselves on the baryon loading 
problem in this paper. In a companion paper, we will 
study in detail how the BH central engine model may 
interpret the phenomenology of GRB prompt emission 
and X-ray afterglow. 
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